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^ ■ Abstract 

I We describe the strong coupling limit {g oo) for the Yang- 

\Q • Mills type matrix models. In this limit the dynamics of the model is 

. reduced to one of the diagonal components which is characterized by 

I a linearly confining potential. We also shortly discuss the case of the 

' pure Yang-Mills model in more than one dimension. 

p 

1 Introduction 

D I The development of string and gauge theories is characterized by their strong 

'-^ I inter-relations. The most intriguing result of this interaction is, probably, the 

AdS/CFT conjecture [T1I2] (see jS] for a classical review of the subject). This 
^ ■ conjecture relates the string theory in the Anti-de Sitter background on the 

one side with the conformal theory on the Minkowski space-time on the other. 
The Minkowski space-time of the conformal theory in this case is related to 
the (conformal) boundary of the Anti-de Sitter space. 

This conjecture relates a weak coupled model to a strong coupled one and 
vice-versa, which is a true Ising-type duality. Once proved, it would have an 
immense predictive force, e.g. for describing the strong coupled dynamics of 
both strings and gauge fields. On the other hand, it is clear, that for a direct 
proof, one needs to know the strong coupled behavior of at least one of these 
models (in addition to the weak coupled one for the both). A considerable 
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progress was achieved in recent years on the way of indirect proofs of the 
correspondence (for a recent review see e.g. 

On the other hand, in spite of difficulties in the description, it seems, that 
the strong coupled regime of the gauge models is the most natural regime 
realized in the Nature at the most common (i.e. low) energies. Perhaps, 
the most success in the description of the strong coupled gauge theories 
was achieved in the framework of the lattice formulation.^ An important 
problem of this approach, however, is that the continuum limit of strong 
coupled systems is problematic and it is difficult to separate the real physical 
effects of the strong coupling from the artifacts of the lattice description. 
Therefore, it would be important to have a strong coupling approach not 
related to the lattice discretization. In the present work we attempt to move 
into this direction. 

Although at the end we also consider the Yang-Mills model, the main 
subject of this paper is the BFSS type matrix model alias Yang-Mills me- 
chanics. Yang-Mills type matrix models appear in both contexts of string and 
gauge theory. Thus, BFSS and IKKT [7j matrix models were proposed 
to describe, respectively the "zero"- and "minus-one" -brane configurations 
in the nonperturbative string approach (M-theory). They can be obtained 
as dimensional reductions to, respectively, 1+0 and dimensions of the ten 
dimensional super Yang-Mills model. (See e.g. [Hj for a review.) 

The plan of the paper is as follows. In the next section we shortly intro- 
duce the matrix model. Then we consider the g —>■ oo limit of the matrix 
model. First, as a warmup we consider what we call a strong limit. In this 
limit we do not consider the contribution from the high frequency modes. It 
leads to a model for the diagonal components where all fields are statistically 
confined (condensed) to a single value. Next, we consider a more refined 
weak limit where we take into consideration the above higher modes. This 
leads to a dynamically nontrivial model for the diagonal components which 
are interacting by linear attracting potential. In addition this model appears 
to be semi-classical as g goes to infinity. At the end of this section we discuss 
a possibility for a systematic expansion at large coupling. 

At the end we discuss the possibility to extend the analysis to the Yang- 
Mills model. 

-'^A good reference for the lattice approach to gauge theories is given by [5]. 
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2 Matrix model 

Consider the matrix model (Yang-Mills mechanics) which is described by the 
following classical action: 

S = J dttr{i(VoX„f + ^[X„X,]2}, (2.1) 

where Xa are D, a = 1, . . . D time dependent Hermitian N x N matrices 
while g is the gauge coupling. The covariant time derivative is defined by 
the use of the (non-dynamical) temporal gauge field A = Aq, 

VoXa = Xa+[A,Xa]. (2.2) 

The role of the gauge field is to impose the Gauss law constraint [Xa, VoXa] = 
which provides the gauge invariance of the action with respect to the time- 
dependent U(A^) gauge transformations, 

Xa^U-\t)XaU{t), A^U-\t)AU{t) + U-^U (2.3) 

where U{t) ElJ{N). Other features of the model include: 

• Invariance with the respect to shifts by a constant scalar matrix 

Xa^Xa + Ca- I (2.4) 

Restricting the gauge group to SU(A^) removes this degree of freedom 

• Invariance with respect to the (target space) rotations, 

Xa ^ Aa%, (2.5) 

A gSO(D) 

• In the case of D = 10 eq. 1)2.11) represents the bosonic part of the 
supersymmetric BFSS matrix model 0, 

j dttr|i(VoXj2 + ^[x„X,]2 + V;VoZ^ + ^r"[X„V^]}, (2.6) 

where ip is the fermionic N x N matrix with 10 dimensional Majorana- 
Weyl fermionic indices. 

For the matrix model under consideration one can formulate a pertur- 
bative expansion in terms of the powers of the gauge coupling g similar to 
the perturbative expansion of the Yang-Mills theory. In what follows we will 
not discuss this type of perturbative expansion but refer the reader to the 
appropriate Yang-Mills perturbation theory literature instead. 
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3 Spontaneous symmetry breaking at strong 
coupling 

It is expected that the strong couphng hmit g oo imphes the commuta- 
tivity of the matrices Xa, 

[Xa,X,]=0. (3.1) 

Indeed, as g goes to infinity the path integral contribution of configurations 
with non-zero commutator are exponentially suppressed. 

Since this is the case, one can diagonalize simultaneously all the matrices 
Xa, whose eigenvalues would then correspond to the coordinates of some 
branes. In this case one can say that in the strong coupling limit the branes 
can be localized. (Beyond this limit they are fuzzed by the strings by which 
the branes interact.) 

Let us consider the above (7 — > oo limit in more details. For this let us 
split the matrix degrees of freedom X^^ into the diagonal part: 

Xa = diagXa, (3.2) 

and the remaining off-diagonal one: 

Za = Xa- Xa- (3.3) 

This splitting seems somehow abusive, since it does not respect the gauge 
invariance ()2.3|) . In fact, it corresponds to a particular choice of commutative 
background among gauge equivalent ones. This choice breaks spontaneously 
the U(A^) symmetry down to U(l)^. At the same time Za can be treated as 
a perturbation above this background. 

The spontaneous breaking of the symmetry is always associated with 
the zero modes corresponding to different gauge equivalent choices of the 
background^. An appropiate choice of SU(A^) gauge apparently solves this 
problem since it restricts the allowed perturbations of the vacuum to the 
transversal direction. The unbroken gauge symmetry as well as the possibil- 
ity to fix the gauge depends strongly on whether the diagonal background is 
degenerate or no. Although the exceptional configurations with the degener- 
ate background may in principle contribute (and even dominate) in spite of 
zero measure, we so far neglect this issue and consider in rest of this paper 
the general position point: where all diagonal eigenvalues different (as 

D-dimensional vectors) . 

^In the present case this is the symmetry: Xa U^^XaU. 
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4 g ^ oo: the strong limit 

On can define different strong coupling limits depending on the relation of 
the coupling with other parameters (like N or the cut-off). In this section 
we consider the strong limit: This limit assumes that the model is UV- 
regularized and the limit (7 ^ cxd is taken prior to removing the regularization. 
Technically, this means that one can drop in this limit the time derivatives 
if they come with a factor vanishing in the limit g ^ 00. In contrast to 
this, the weak limit which is taken after the removal of (or eventually not 
imposing at all) the regularization is discussed in the next section. 

In the non-degenerate case the whole U(A^) gauge group is broken by the 
diagonal component of the background down to U(l)^. The infinitesimal 
gauge transformation of the background is given by 6xa = and 6za = 
[xa, u] + [za, u]. This is very similar to the ordinary gauge transformation in 
the nonabelian Yang-Mills theory if the role of the partial derivative operator 
da is attributed to the commutator [xa, ■]■ In the complete analogy with this 
one can fix the gauge by imposing the Lorenz gauge condition^: 

F^,, = [Xa, Za] = 0. (4.1) 

The Faddeev-Popov determinant corresponding to the gauge fixing condition 
()4.1|1 is given by 

1 

(4.2) 

where the prime denote that the product extends over the distinct indices m 
and n only. Formally, the determinant is different from zero (which is impor- 
tant for the implementation of the gauge condition) when all x-eigenvalues 
are given by distinct points x^, n = 1, . . . N. 

All above can be appropriately formalized in the quantum theory by 
adding the gauge fixing term and the Faddeeev-Popov determinant in the 
(Euclidean) partition function which takes the form 

Z = j [dx][dz][dA]At\x)exp{-{S+'4tT[xa,Za]')}, (4.3) 

where we used so called "alpha-gauge" (with a = g^) implementation of the 
gauge fixing rather than the "delta-function implementation". Note, that 

■^Admissible gauge fixing and corresponding Faddeev-Popov determinants are discussed 
in the classical book on gauge theories 0. 



x)=n 



time 



n' 
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since the introduction of the gauge fixing condition ()4.1|) one cannot anymore 
impose any further restriction^ on the gauge field A which should remain in 
the action. 

Now we are ready to take the limit g ^ oo and separate the leading 
contribution in this limit. There are several ways to do this, which, naturally, 
lead to the same result. Let us consider the following one. Let us substitute 
the variables Za by the rescaled ones as follows 

Za gZa- (4.4) 

Then, the matrix action takes the following form: 



S=- j dt {^\Xa^ + ^Za^ + Za]Xa + \[A, (x, + ^^a)]' + 

{Xa + -gZa)]Za + \{[Xa, ^fe] " [Xfc, Za] + ^[^a, ^b])^ + Zaf]^ . (4.5) 

As we are taking the strong g oo limit, we should discard all terms formally 
vanishing in this limit. Thus, the leading part of the action becomes 

Sg^oo = - j {¥-^ + |[Aa:a]' + \[xa.zi,f) . (4.6) 

The action is quadratic in the gauge field A as well as in the off-diagonal 
field Za- Integrating in both A and Za, one gets the factor coinciding with 
the Faddeev-Popov determinant at the power —{D + l)/2. The partition 
function then reads, 



da;A(°°)(x)-(^-i)/2 



exp <^ dt^x'}, (4.7) 



which appart from the determinant factor in the measure corresponds to a 
free particle partition function. 

The modification of the measure in ()4.7|) signals the confining of the 
eigenvalues Xn to a common value which itself is a subject to free motion. 
Indeed, in the case of only two eigenvalues the path integral ()4.7|1 reduces to 
(see the Appendix 



j [d^F] [d^yy-2(^-i)]e'^"<^'"+^^'), (4.^ 



where y is the distance between branes while Y is the free moving "center 
of mass" . Consider the |/-measure locally at the instant t: d^y{y^)-^^\t) = 

^Except for the vanishing of the diagonal part of A, Ann = 0. 
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dQ£)drr~^^^ . Integration with such a measure is divergent at r = ^^/y^ = 
unless the integrand vanishes quickly enough as y approaches the origin, 
which is not the case for slow y modes. Statistically this means that con- 
figurations with small y"^ produces a contribution to the partition function 
which is infinitely larger than the contribution of all the configurations with 
larger values of y"^. Therefore, under the normalization the configurations 
with nonzero y"^ will get zero expectation values. One can see also that the 
conclusion is very sensitive to the power of A^°°\ Thus, if the power were 
e.g. —{D — l)/4 no such statistical confinement would occur. 

It may appear however that this simple estimation of (7 — > cxd is too rough 
and one must weaken the limit allowing the contribution of higher frequency 
modes. We come to this in the next section. 

5 ^ ^ 00: the weak limit 

Consider the stationary points of the action ()2.1|) i.e. the solutions to the 
equations of motion. There is a class of static solutions to the equations of 
motion given by constant commuting matrices Xa- We can assume that these 
matrices depend adiabatically on time. One can consider perturbations about 
this background. The perturbation is given by the off-diagonal part Za as well 
as by the fast diagonal modes. The diagonal modes do not contribute at the 
one-loop level since there are no nonlinear terms in the action corresponding 
to diagonal-diagonal interaction. As a sequence, we can neglect the fast 
diagonal fiuctuations and consider only the adiabatic modes. 

Therefore, consider the contribution of the off-diagonal modes as well as 
of the auxiliary (gauge and ghost) fields and evaluate their contribution in 
the one-loop approximation in 1/g expansion. Throughout this section we 
use the Euclideanized version of the theory. 

To proceed with the evaluation let us fix the gauge by adding the following 
gauge fixing term to the Lagrangian: 

L,,f,=-tT(^^A' + l[Xa,Z,]'y (5.1) 

The variation of the gauge fixing condition gives the Faddeev-Popov opera- 
tor, 

MppU = doVoU + [Xa, [{Xa + Za),u]], (5.2) 

whose determinant App = det Mpp is the Faddeev-Popov determinant which 
we have to use together with the condition (jOJ). In the one loop approx- 
imation no contribution will come from A- and z-dependent terms in the 
Faddeev-Popov operator. Therefore, in what follows we will discard these 
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terms. As a result, the Faddeev-Popov determinant restricted to one loop 
relevant terms takes the following form 



Fp|(i loop) = n + = ^2ix), (5.3) 



A, 



where = (x^ — x^)^ is the square distance between n-th and m-th branes 
and the prime denotes that the product is taken for distinct m and n. 

Let us turn to the action ()2.ip . The matrix model action can be rewritten 
in the form as follows, 

•S* = — ^ dt ^2 (^n) ~^ 2g^\^mn\ ~^ 2g^\^mn\ ~^ '^Cmn.Cmn 

2^mni.\^mn\ ~^ l^mnl ~l~ ^mn^rmi) "!"•••) (^'4) 

where the dots stand for the terms not contributing at the one loop level 
(e.g. terms which are higher than the second order in A and z). 

After the integration over the gauge field A, the off-diagonal component 
z and the ghosts c and c the partition function takes the form. 



D-l 



[dx]A2 2 (x)e^''*i"' (5.5) 

As it can be seen, the problem is reduced to the computation of the de- 
terminant A2, of an elliptic differential operator. Let us use the (^-function 
approach to do such a computation^ (see e.g. P^). According to this ap- 
proach, the logarithm of the determinant of an elliptic operator D is given 
by the (minus) derivative of the ^-function, 

IndetD = -C^(0), (5.6) 

where the function Cd{s) is defined as the analytic continuation of the series, 

1 1 /"°° 

&W = EF = fwi dpp-tre-- (5.7) 

A 

The trace tre~^^ can be written as 

tre~''^ = / dtKD{t,t]p), (5.8) 



similar computation for constant diagonal modes was done in and its phe- 
nomenological implications were explored in (15) . I thank Amir H. FathoUahi for pointing 
my attention to these papers. 
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where K£,{t' ,t"] p) is the Heat Kernel for the operator D. It is the solution 
to the Heat Equation 

d,K{t,to;p) = -DK{t,to;p), (5.9) 
with the initial conditions given by 

K{t,to;0)=6it-to). (5.10) 

In the case at hand D = {^^d^ + r^^j and the solution for the Heat 
Kernel is given by 

Kit ,t ]p) = ^=e^]i[ r„nP • (5-11) 



V 4p 

Since the time integral in the r.h.s of the equation ()5.8|1 diverges for t G 
(— oo, +oo) it is useful to put the system in the time box interval r. Beyond 
its regularization function the r plays another important role, namely, that 
of being also the adiahaticity box. Roughly speaking, the r-interval is the 
"dt" for the adiabatic time "t" . 

The C-function for the time interval r is then given by 







mn 



Computing the derivative of ()5.12j] and taking the limit s ^ we obtain: 

-CoiO) = grY! V^^. (5.13) 

mn 

Summing over the all adiabatic boxes we get: 

A2~(x) =e-^^-^<^*^™v^, (5.14) 

where we can even drop the prime from the sum. 

Therefore, the low energy effective action for x'^ takes the form. 



\ mn 



Sg^^ = dt \ -\xi - \g{D - 1) > ^ V(x„ - x„)2 . (5.15) 
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As one can see, the action ()5.15|) corresponds to a system with strong 
hnear confinement of the particles. In spite of its terrifying appearance the 
hmit g ^ oo corresponds to nothing else then the semi-classical limit. Indeed, 
passing to a rescaled x^, 

(5.16) 

transforms the partition function (jS.lSj) to the following semi-classical form 
Z = J [dx]e'"^''-'^" ^ ^ (5.17) 

where plays the role of inverse Planck constant h^^. In fact, the above 
rescaling introduces a renormalization of the brane coordinate. Its meaning 
is that the nontrivial dynamics corresponds to large (in the old scale) brane 
separations. Therefore, the natural scale of the brane dynamics is given in 
terms of the attraction force (tension) acting on the branes. 

A remark on the systematic expansion 

A trick can be used to modify the value of the coupling constant (and even 
to invert it). 

We can consider the model at the finite temperature T = 1//5. The finite 
temperature implies that the action in the path integral is computed for the 
Euclidean time interval < t < (3 with periodical boundary condition for 
the fields. A simple dimensional analysis that the following rescaling 

P^P/X', g'^g'X', (5.18) 

changes the partition function by a constant multiplicative factor only, which 
can be absorbed in the measure. Indeed, making the substitution X \X 
one gets ()5.18j) . Now taking A arbitrarily small one can make g small as well, 
e.g. equal to g~^. At the same time f3 goes to infinity i.e. the theory rolls 
down to zero temperature. 

Unfortunately, because of different scaling properties, this trick can not 
be used in the case of Yang-Mills theory in more than two dimensions. 

6 The Yang— Mills model 

It is tempting to apply the above analysis to the SU(A^) Yang-Mills model. 
Let us enumerate the modifications that occur when passing to the pure 
Z)-dimensional Yang-Mills model: 
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Instead of the determinant ()5.3|) one should compute the determinant 
of the D-dimensional differential operator 

D = j,dl-rl^, (6.1) 

where the diagonal modes are described by the Abelian gauge fields 
'^^mn = ('^/T' ~ '^^Y ■ Also since the gauge group is SU(A^) the 
center of mass is fixed: 

Y.al = Q. (6.2) 



Heat Kernel: 



• The one loop contribution is given by:^ 

Leff = \{F;,f -g''iD-2)J2 VUa), (6.4) 

mn 

where FJ^^ = d^a^ — d^a^ and 

™" " (4-)-/^ \^ (logrL - h{D/2)) , D-even ^'"'^ 

where h{k) is the fc-th harmonic number: h{k) = Yl'i=i ^/^- (Note also 
that the F-function is regular at negative half- integer points.) 

As it could be seen, for D > 2 {D = 2 is dynamically trivial) one can 
rescale the fields 

D 



a;^gD-2a;, (6.6) 

2D 



and get a common factor g'^-^ in front of the effective action. For D > 2 
this factor vanishes in the limit g ^ oo, which means that in this case 
the quantum fluctuations are strong. As one can note, the the qulitatative 
behaviour of the effective models depends on dimension. Thus, in dimensions 
D = Ak and 4fc + 1 for non-negative integer k, the strong attractive force 
binds al a" together, while for D = 4k + 2 and 4/c -|- 3 the repulsive force 
keeps them appart at infinity. The common feature is that in this situation 



^The factor {D — 2) instead of {D — 1) as in the previous sections appears because the 
gauge field Aq is now counted as one of the fields. 



11 



we are not able to catch any nontrivial dynamics beyond the fact that all 
diagonal values are confined to zero or infinity. 

A much more serious problem is that ioi D > 2 the higher loop contri- 
bution is not suppressed at large g unless an UV-cutoff (A < oo) is used. A 
nontrivial contribution can be then catched taking the double scaling limit 
with g ^ oo and A ^ cx3. A more detailed analysis would give the answer 
whether this is possible. 

7 Discussion 

In this paper we considered the strong coupling limit of the matrix model. 
It is shown that the modes which survive in this limit are described by 
a system of linearly interacting particles. As coupling goes to infinity the 
system becomes semi-classical g^ playing the role of inverse the inverse Planck 
constant h~^. The scale at which the dynamics takes the semi-classical form 
is given then by the string tension or the coefficient of linear interaction. The 
analysis is performed at the one-loop level. It seems rather possible that a 
systematic expansion in the inverse powers of the coupling constant can be 
constructed in addition to the standard small coupling constant expansion. 

It is also very tempting to apply the 1/ (^-expansion to the Yang-Mills 
theory. The one-loop technique can be easily extended to the ordinary Yang- 
Mills model. In the case of the dimensionality higher than two the diagonal 
component is not anymore semi-classical and most probably does not decou- 
ple. The implications of this are not yet clear. There are however, resources 
we did not use which are given by the large and UV cut-off scaling. Tak- 
ing a correlated limit of large g, N and A one may hope to get a non-trivial 
content for the expansion, e.g. by tuning the background. 

Another important issue we left beyond our consideration regards the 
exceptional configurations with some r^n = 0. As the effective parameter of 
the expansion is 1/grmn the expansion fails if some r^n ^ g~^- Important 
point is the statistical weights of such configurations. An estimate can be 
done by the computation of the average separation f. When the average 
separation is nonzero f > 0, it is clear, that one can trust the approach. In 
the case of pure Yang-Mills model, however, it seems that it either vanishes 
or is infinite according to the dimension. 

In the case of branes at close distance the 2x2 matrix block corresponding 
to respective eigenvalues is not decoupled and one should consider the entire 
matrix dynamics similarly to what is done in the non-commutative gauge 
theory [HI [l2l HHj . As it was found this dynamics is a stochastic one. 
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A Example: The tale of two branes 

Consider the case of two branes. In this case the action can be written in 
the following form: 

+ [a{y ■ z — z ■ y) — a{y ■ z — z ■ y)] — {ai — 02) {z ■ z — z ■ z) 
+ \p2 {a\ — 02) [ay ■ z + y ■ za) — 2aa{z ■ z + y^) + z'^a'^ + a^z^ 

— \z-ziax — a'lf' 
- {2y^z ■z+{z-zf - z'^z^) - 2y^{cc + c*c*), (A.l) 

where 2x2 matrix Xa is given by the following component structure, 

while the gauge field matrix A is given by the components: 

A=f^i (A.3) 

and two complex conjugate components of the ghost-anti-ghost are used. All 
diagonal components are real while the off diagonal elements are complex. 
The dot in (jA.ljl indicates the inner product with respect to the index a = 
1,...,D. 

The first four lines of ()A.1|) are the contribution from the kinetic term 
while the third line comes from the commutator term together with the gauge 
fixing term and Faddeev-Popov determinant for the gauge [x, z] = 0. 

Let us make the following substitution: 

z gz. (A.4) 

After the rescaling one can split the Lagrangian ()A.1|) in the leading term 
and perturbation in 1/g. The leading part of the Lagrangian looks as follows, 

Cg^oo = \y'^ + \y^ -2y'^{aa + z-z + cc + c*c*). (A.5) 

13 



All fields with the exception of the y and the free Y become non-dynamical 
in the limit g oo and the Lagrangian ()A.5|) is quadratic this fields. There- 
fore, integration of and a leads'' to the partition function of the 
form (gUj), 

Z= f J]d^Frf^i/b2(t)]-^+ie'^'<^^'+^*'). (A.6) 

t 

As in the case of ()4.7|1 the measure in eq. ()A.6|1 is singular as ^ 0. 
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